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(N ! ABSTRACT 

m 

Context. We have observed again two stars inside the ZZ Ceti instability strip that were previously classified as not- 
observed-to-vary (NOV) by Mukadam et al. (2004) and found them to be low-amplitude variables. Some evidence points 
to a pure ZZ Ceti instability strip; other evidence contests it. 

Aims. The two stars previously classified as NOV have Sloan Digital Sky Survey (SDSS) spectroscopic effective tem- 
peratures that place them inside the ZZ Ceti instability strip, and they were "contaminating" the strip as constant 
stars, which could indicate that the instability strip was no longer a simple evolutionary stage. A pure instability strip 
indicates that pulsation is a normal phase which all DAs must go through. 

Methods. We used effective temperatures derived from SDSS optical spectra by comparing them with model atmo- 
spheres to look for pulsators through time-resolved photometry and stars previously classified as NOV. 
Results. Our new results indicate, but do not prove, a pure instability strip, because there are still other NOV stars 
that need to be observed again. Additionally, we have discovered five other ZZ Ceti stars based on their effective 
C$ • temperatures. 

Key words. (Stars): white dwarfs, Stars: variables: general, Stars: oscillations 

& ' 1. Introduction the DA white dwarf stars with temperatures placing them 

_,, . . . . .... . , . . inside the instability strip do not show pulsation above their 

Ihe observational DA instability strip (ZZ Cetis ranges , , ,. ,. ., n , , , , , 

„ ; . ,—. , J . ,i„.„ \ nn r^-rr detection limit, so were called not-observed-to-vary (NOV), 

m effective temperature (T e g) from 12 270 to 10 850 K A/r ,. ' • , , , u 

lr ., . _, 1 -A . nnnr - t»t t i i Many questions were raised based on these new observa- 

Giannmas, Bergeron, & lontame 2005; Mukadam et al. T .v • . i.,-, . • 9 ,,ri. i .v _i ■ , 

X__ .% ' , , i • , , , /t-, • , i tions. Is the instability strip purer Which are the physical 

2004). rhe current best theoretical models (lontame et al. u ■ j-v, + i. i„ ±- 9 ™tu j 4-r,- 

t-, 11 a m 1 n n-i i , mechanisms that can prevent pulsation/ What does this 

2001: Bradley 1996; Arras Tbwnsley, & Bildsten 2006 pre- particular gample haye that ig different from the old Qne? 

diet that the ZZ Ceti instability strip should be pure: all the Can the H j be extremel tMn so that pulsations win 

DA white dwart stars, with temperatures within the msta- , , o ri 1 ±- c u 1 

.... , •,, ,1 t 1 not happen r Can even weak magnetic helds prevent pulsa- 

bihty strip hmits should pulsate, with a small dependency ^ Dq thfi N0Vg haye unfavorable viewin an les? 

on mass. Ihe theory agreed extremely well with observa- 
tions until the search for new pulsators was extended to If the instability strip is really contaminated with non 
white dwarfs discovered by the Sloan Digital Sky Survey Pulsators, the models for DA white dwarfs and the results 
(SDSS), which reaches a fainter population (16 < g < 20) derived from these models need to be revised. Before taking 
than the old sample. Using accurate temperature determi- such a crucial but drastic step, it is important to make 
nation either from optical or ultraviolet spectra, the insta- sure the observational results are correct. There are several 
bility strip using only the old bright sample did not con- wa y s to falsif y a contaminated instability strip. First, the 
tain any non pulsator (e.g. Bergeron et al. 2004; Gianninas, temperatures of the SDSS sample are derived from spectra 
Bergeron & Fontaine 2005) with a signal-to-noise rate (SNR) <35. In this scenario, 
However, in the searches' conducted by Mukadam et al. the temperatures and masses of the non variables could 
(2004) and Mullally et al. (2005), a substantial fraction of Slm * be inaccurate and the stars actually he outside the 

instability strip (Gianninas, Bergeron, & Fontaine 2005). 
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these uncertainties, and the searches should be extended to 
a wider range of temperatures (e.g. the odd pulsator WDJ 
2350-0054 at 10 350 K in Mukadam et al. 2004). Second, 
pulsating white dwarf stars can appear to be extremely low- 
amplitude pulsators for a few hours due to the beating of 
many excited modes (Kleinman et al. 1998). A third but not 
the least important possibility is that variations could have 
been, in fact, smaller than the published detection limits, 
which vary from 2 to 9mma (eg. Kanaan et al. 1992). 

The only way to solve the question about the purity of 
the ZZ Ceti instability strip is to observe these stars again 
through higher SNR spectroscopy and time-resolved pho- 
tometry. The higher SNR spectra allow a more accurate 
determination of T e g and log g (Kepler et al. 2006). In this 
paper, we present our first results on the re-observations 
of the time-resolved photometry of two stars classified as 
NOVs, finding variability in their light curves. As every 
NOV we observed is in fact a low-amplitude variable, we 
show strong observational evidence of a pure ZZ Ceti insta- 
bility strip. 

2. Observations 

We observed our target list with the SOAR 4.1-m tele- 
scope, using the SOAR Optical Imager, a mosaic of two 
EEV 2048x4096 CCDs, thinned and back-illuminated, with 
an efficiency around 73% at 4000 A, at the Nysmith focus. 
We observed in fast readout mode with the CCDs binned 
4x4, to decrease the readout -|-write time to 6.4s and still 
achieve a 0.354" /pixel resolution. The exposure times were 
30 s, and we observed up to 4 hours. All observations were 
obtained with a Johnson B filter to maximize the amplitude 
and minimize the red fringing. 

We also observed with the Otto Struve 2.1-m telescope 
at McDonald Observatory, using the Argos camera. This 
camera is mounted at the prime focus of the telescope, it has 
a frame transfer CCD, and uses a red cutoff BG40 filter to 
reduce the scattered light from the sky (Nather & Mukadam 
2004) . The integration times were 5 s and the observations 
spam for up to 1.4 hours. 

Extending our search for new pulsators, we used the 
1.6-m telescope at Observatorio do Pico dos Dias, LNA, 
in Brazil. For these observations, we used a frame transfer 
CCD, with a quantum efficiency of 60% around 4000 A. We 
observed with no filters because the CCD used has no de- 
tectable fringing and the g-mode pulsations in ZZ Ceti stars 
are coherent at all optical wavelengths (Robinson, Kepler, 
& Nather 1982). The integration times were from 40 to 50 s 
and spanned up 2.5 hours. In Table [T] we show the journal 
of observations. 

For each telescope, we chose the filter set up for which 
the SNR were optimal. 

The technique used to detect variability is differential 
time series photometry, comparing the targets with the 
other stars in the same field to minimize the effects of sky 
and transparency fluctuations. For each run, we extracted 
light curves using the hsp scripts developed by Antonio 
Kanaan for IRAF, with weighted aperture photometry, us- 
ing different aperture sizes and selecting the light curve 
with the highest SNR in the Fourier transform. 

For the stars previously classified as NOVs, we had to 
observe them for a minimum of three hours. This mini- 
mum time was chosen to avoid destructive beating between 
modes (if there are any that are excited in the star) and 



to be able to reach a lower amplitude limit than previous 
observations. The NOVs we observed had detection limits 
above 3-4 mma. 

In order to have some objective criterion for determining 
which peaks are real in the discrete Fourier transform, we 
adopted an amplitude limit such that a peak exceeding this 
limit only has a 1/1000 probability of being due to noise 
(false alarm probability). For these continuous data sets, 
peaks above 3 (A) (3 times the square root of the average 
power) have a probability smaller than 1 in 1000 of being 
noise. We checked the detection limit for each star by Monte 
Carlo simulation, randomizing the acquisition times. 

To know if another peak in the Fourier transform is an 
intrinsic periodicity of the star or whether it is only due to 
the spectral window, we subtracted from the light curve, 
the sine curve with the same amplitude, period, and phase 
information as the peak selected in the Fourier transform. 
After subtracting it from the light curve, we recalculated 
the Fourier transform to verify whether the sine curve was 
correct. Then, we repeated the procedure for the remaining 
periodicities. 



3. New variables 

The stars SDSS J030325. 22-080834.9 and SDSS 
J085325.55+000514.2 were listed in Mukadam et al. 
(2004) as NOV objects with detection limits of 3-4 mma. 
In Fig. [T] we show the Fourier transform of the new 
pulsators previously classified as constant stars, both 
being low- amplitude pulsators, with amplitudes below the 
previous detection limit. 

We show in Fig. [5] the Fourier transforms of the to- 
tal light curves of the new stars discovered to be pulsators 
that have never been searched before. These stars were in- 
corporated in our candidate target list based on the same 
criterion: their T e g, derived from fitting the SDSS optical 
spectra to Koester's synthetic spectra derived from model 
atmospheres. We used a model atmosphere grid similar to 
the one described by Finley, Koester, & Basri (1997), but 
larger and denser. Within the T e g range from 11000 to 
12 000 K, we have a 90% probability of finding pulsators 
(e.g. Mukadam et al. 2004). 

Table [2] shows a summary of the physical properties of 
the new pulsators derived by fitting the optical spectra up 
to 7000 A to Koester's models in the same way as described 
in Kleinman et al. (2004) and Eisenstein et al. (2006). We 
measured the x 2 on the difference between the observed 
spectra and the models, allowing for a reflux of the models 
according to a low-order polynomial. This step was made 
to correct for effects of unknown reddening. We used the 
first seven Chebychev polynomials in a linear wavelength as 
our refluxing basis. The uncertainties quoted in the table 
were estimated as the internal uncertainty of the fit only, 
not considering that T c g and logg are correlated. These 
estimates also do not consider the changes in T e g during 
the pulsation (see Sect. 4.1 for a detailed discussion) or 
external uncertainties, which are dominant. 

We analyzed the light curves of each one of the new 
pulsators to identify their periodicities, necessary in fu- 
ture seismological analysis for studying their properties and 
structure. The list of detected periods in these stars is 
shown in Table [3l 
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Table 1. Journal of observations with 2.1-m telescope at McDonald Observatory, 4.1-m at SOAR, and 1.6-m at OPD. 



Star 


Run start (UT) 




AT (hr) 


# points 


Telescope 


SDSS J000006.75-004654.0 


2005-08-04 07:29 


45 


2.5 


199 


1.6-m 




2005-08-06 05:23 


45 


1.6 


127 


1.6-m 




2005-08-07 05:53 


50 


1.2 


88 


1.6-m 


SDSS J030325.22-080834.9 


2005-12-05 02:59 


30 


3.0 


301 


4.1-m 




2005-12-06 02:29 


30 


3.5 


347 


4.1-m 




2005-12-07 01:10 


30 


4.0 


400 


4.1-m 


oJJDD JU0OOZO.OO + UUUOI4.Z 


zUUo-lz-Uo Uo.lU 


oU 


6.6 


zzy 


/11m 

4.1-m 




2005-12-09 09:50 


30 


3.5 


170 


4.1-m 


SDSS J085507.29+063540.9 


2005-12-11 05:02 


30 


2.8 


276 


4.1-m 


SDSS J091635.07+385546.2 


2006-01-07 06:35 


5 


1.3 


921 


2.1-m 




2006-01-08 08:46 


5 


1.4 


969 


2.1-m 


SDSS J165020.53+301021.2 


2005-08-02 00:37 


40 


2.2 


191 


1.6-m 


SDSS J233458.71+010303.1 


2005-08-06 06:14 


50 


2.2 


306 


1.6-m 




2005-08-12 00:36 


30 


1.5 


198 


4.1-m 



Table Notes. AT is the total length of the each observing run and t cxp is the integration time of each exposure. 



Table 2. Physical parameters derived by comparing the whole SDSS optical spectra with a grid of model atmospheres 
and SDSS magnitudes. 



Star 


T efi (K) 


logs' 


g (mag) 


Mass (Mq) 


Tc« (K) «P» 


WDJ 0000-0046 


10 880+110 


8.32+0.09 


18.84 


0.81+0.06 


11310 


WDJ 0303-0808 1 


11400+110 


8.49+0.06 


18.74 


0.92+0.04 


10 960 


WDJ 0853+0005 


11750+110 


8.11+0.06 


18.23 


0.68+0.04 


11710 


WDJ 0855+0635 


11050+50 


8.43+0.03 


17.25 


0.88+0.02 


11 140 


WDJ 0916+3855 


11410+50 


8.10+0.03 


16.56 


0.67+0.02 


11410 


WDJ 1650+3010 


11 100+90 


8.76+0.08 


18.11 


1.07+0.04 


11620 


WDJ 2334+0103 


11400+210 


7.99+0.14 


19.24 


0.60+0.08 


10 930 



Table Notes. We derived the masses using T c h and logg obtained from optical spectra and comparing these values with C/O white 
dwarf evolutionary models (Wood 1995). We also calculated the T e g using its relationship with weighted mean pulsation period 
(see Table|3]) derived by Mukadam et al. (2006). The quoted uncertainties are the internal uncertainties due to fitting procedures 
only. 1 Values derived from Gemini spectra are T eff = ll 960+160 and log g=8. 305+0.017 using LPT fit and T cff =11 420+110 
and log g=7.821±0.016 when the fit procedure included the whole spectrum (Kepler et al. 2006). 



Table 3. Periodicities detected in our discovery runs for 
the new ZZ Ceti stars. 



Star 


Period (s) 


Amplitude (mma) 


WDJ 0000-0046 
WDJ 0303-0808 
WDJ 0853+0005 
WDJ 0855+0635 
WDJ 0916+6108 

WDJ 1650+3010 
WDJ 2334+0103 


611.42; 584.84; 601.35 
707; 1128 
264.35 
850; 433 
485.09; 447.70; 
238.10; 747.20 
339.06 
923.15 


23.00; 15.92; 8.97 
4.1; 3.5 
3.99 
44; 15 
32.89; 14.44; 
10.76; 9.06 
14.71 
40.37 



4. Discussions 

4.1. Purity of the instability strip 

In Fig. [3] we show the updated SDSS ZZ Ceti instabil- 
ity strip, including the NOVs reported by Mukadam et al. 
(2004) and Mullally et al. (2005) that still need to be re- 
observed. 

At this point it is necessary to discuss the star WDJ 
2350-0054, which has a published T cS = 10 350 ± 60 
(Mukadam et al. 2004) , indicated by an aster- ix in Fig. [3] 
It is a pulsating star outside the instability strip with T e s 
significantly cooler than the red edge. Because there is no 
reasonable explanation for such a low temperature, we will 
not consider this star as the red edge of the instability strip. 
There are three independent SDSS spectra for this star, 



all of them providing a lower T e g than does the red edge 
when the whole spectrum is used in the fitting procedure. 
However, our line profile technique (LPT) fit, which has a 
larger scatter in the T e g determination for the SDSS spec- 
tra because of the lower SNR in the lines, the resulting T e s 
brings this particular star inside the instability strip. High 
SNR spectra are necessary to reliably determine its physical 
parameters. 

We detected variability in two stars previously classified 
as NOVs. There are still 19 NOVs inside the instability 
strip, which need to be observed again. A strong argument 
for a pure instability strip is that every star we re-observed 
for variability turned out to be a pulsator. 

Kepler et al. (2005) discuss the external uncertainty 
measured from SDSS duplicate spectra, showing the un- 
certainty is around 300 K in T e g and 0.21 dex in log g, 
a significant portion of the size in T e g of the instability 
strip (~ 1400 if). There is also another important effect 
on the contamination of a pure instability strip with non- 
pulsators: temperature variation during a pulsation cycle. 
For a low-amplitude low I pulsator, the changes should be 
rather small, around 50 K, but for larger amplitude pul- 
sators, temperature variation can reach 500 K (Robinson, 
Kepler, & Nather 1982). 

The Gianninas, Bergeron, & Fontaine (2005) simulation 
of an LPT fitting shows that uncertainties for the SDSS 
spectra are on the order of 450 K for the average SNR^20 
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Fig. 3. Observational instability strip containing only SDSS stars observed for variability. The filled circles represent the 
ZZ Cetis and the open circles the stars for which variation in the light curve was not observed (NOVs) in previous works 
that need to be re-observed. The filled squares are the new ZZ Cetis reported in this paper, and the filled triangles are 
the stars that were previously classified as NOVs, but we detected variability. The asterisk marks a variable star that, 
according to SDSS T e g, is outside the ZZ Ceti instability strip. 



in our sample. These are external uncertainties of the LPT 
method; however, we used the whole spectra fitting. 

Kepler et al. (2006) report results from Gemini spectra 
with SNR around f00 for four stars from the SDSS sam- 
ple. They show that their new T e g were similar to the ones 
obtained with SDSS spectra. The SDSS quoted uncertain- 
ties are underestimated by 60% in T e g and by a factor of 
4 in logg, and the correlation between T e g and log g is an 
important factor. Their simulation indicates that fitting us- 
ing the whole spectra with reliable flux calibration is better 
than only an LPT fit. 

All these uncertainties added together can falsify our in- 
terpretations of the impurity of the instability strip. Higher 
SNR spectra averaging out pulsations (e.g. Bergeron et al. 
2004), as well as longer exposure time-series imaging, are 
required for lowering the uncertainties in temperature and 
the detection limits of pulsations for all stars within the 
boundaries of the ZZ Ceti instability strip. 



4.2. Changes in pulsation spectrum with temperature 

While the star cools down inside the instability strip, its 
pulsation spectrum changes, as do the pulsation ampli- 
tudes. Evolving from the blue (hot) to the red (cold) edge, 
the periods get longer and their amplitudes larger. For ex- 
ample, the hottest known pulsator, G 226-29 (Gianninas, 
Bergeron, & Fontaine 2005), has a dominant mode at 109 s 
and an amplitude of around 3 mma (Kepler et al. 1995) 
in the optical and no other periodicities besides the com- 
ponents of this triplet. On the other hand, the cooler star 
G 29-38 has 20 identified modes (Kleinman et al. 1998), 
and its pulsation spectra change dramatically from run to 
run. In one season, this star showed no detectable pulsation; 
but its dominant mode is normally around 700 s, and the 
pulsation amplitude for a single mode can reach 60 mma. 
All five of the new variables presented in Fig. [3] are long- 
period high-amplitude pulsators, while the previous NOVs 
have low amplitude, independent of their periods. 
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Fig. 1. Fourier transform of all the data for the stars 
previously classified as NOV WDJ 0303-0808 and WDJ 
0853+0005. The dotted lines represent the detection limit 
above 3 (^4), for which peaks have a 1/1000 probability of 
being due to noise. 
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Fig. 2. Fourier transform of our search for new pulsators us- 
ing the telescopes 2.1-m at McDonald Observatory, 1.6-m at 
Observatorio do Pico dos Dias, and the 4.1-m at SOAR. The 
candidates were selected by their T e ff derived from compar- 
ison of optical spectra with model atmospheres. The dotted 
lines represent the detection limit above 3 (A), for which 
peaks have a 1/1000 probability of being due to noise. 



The qualitative explanation for the change in pulsation 
spectra within the ZZ Ceti instability strip is based on the 
behavior of the convection zone, starting at the base of 
the partial ionization layer, which moves inwards as the 
star cools. The convection zone gets bigger, allowing more 
modes to be excited, as the star progresses to the red edge. 
Because the base of the convection zone is deeper, the ther- 
mal time is longer, so the excited periodicities are rather 
long. In addition, as the position of the convection zone 
is more internal for cooler stars, more energy is needed to 
move the upper layers, which leads to higher amplitudes. 
Note that the timescale for crossing the instability strip is 
almost 1 Gyr. 

Even though we understand how pulsations start, be- 
cause of the increase in opacity at the developing partial 
ionization zone, we do not understand how the pulsations 
stop. All evolutionary models predict a much cooler red 
edge than the coolest ZZ Ceti observed, even though Wu 
& Goldreich (1999) propose that pulsation should stop for 
periods around 1 400 s (as observed) in their models of 
convection-dominated DA pulsations. Kanaan, Kepler, & 
Winget (2002) looked for low-amplitude pulsators at the 
red edge, which would represent when the ZZ Cetis stop to 
pulsate, but they found none. Mukadam et al. (2006) found 
seven low amplitude pulsators with long periods, which is 
too large a number to be explained by inclination effects, 
especially when we consider that different ro's cancel at op- 
posite inclination angles. In our sample, WDJ 0303-0808, 
shows long periods and low amplitude. 

Another reason for low amplitudes occurs for high mass 
stars, when their interiors are partially crystallized, like for 
BPM 37093 (Kanaan, Kepler, & Winget 2005), as the pul- 
sations are restricted to the upper layers above the crystal- 
lized core. However, to test this hypothesis for the low am- 
plitude long period star observed, it is necessary to obtain 
spectra with higher SNR than the SDSS spectra, specially 
in the region of the lines H8 and H9, because they are more 
sensitive to gravity, so measuring T c ff and log g accurately. 

Mukadam et al. (2006), following Kanaan, Kepler, & 
Winget (2002), discuss the relationship between weighted 
mean period and spectroscopic temperature, after studying 
the changes in pulsation properties across the ZZ Ceti in- 
stability strip. They found a linear relation between these 
two physical quantities for all ZZ Ceti stars with available 
temperatures. 

Among the new pulsators, the two stars classified as 
NOVs are low amplitude pulsators, but they have com- 
pletely different main periodicities. WDJ 0303-0808's de- 
tected periodicities are 707 and 1128 s, typical of the red 
edge cool pulsator, while WDJ 0853+0005's main mode is 
around 264 s, typical of a blue edge pulsator. Either WDJ 
0303-0808 is a high mass star, where crystallization already 
occurs in its interior or it is an example of a pulsator stop- 
ping to pulsate. In Tabled we listed T e g and logg derived 
from SDSS spectrum. Using these values and the radius 
calculated by Althaus et al. (2005), we derived a stellar 
mass of 0.93 Mq. Kepler et al. (2006) derived a lower stellar 
mass of 0.83 M Q and 0.52 Mq from Gemini spectra fitting 
with LPT and whole spectra, respectively. As crystalliza- 
tion is important within the ZZ Ceti instability strip only 
for stars with masses higher than ~ 1 Mq , when a signif- 
icant portion of the star will be crystallized, this star is 
most likely an example of a star stopping to pulsate. On 
the other hand, the high logg determined from SDSS spec- 
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tra would indicate that the star WDJ 1650+3010 could be 
crystallized; however, its observed amplitude does not in- 
dicate that. The known correlation between T e s and log g 
determination could again explain this issue, as the deter- 
minations were based on low SNR spectra. 

Using the Mukadam et al. (2006) relation, the T e e for 
the star WDJ 0853+0005 should be 11,700 K in agree- 
ment with the one derived from the low SNR spectra 
11,700 K. For WDJ 0303-0808, the pulsation periods in- 
dicate 11,000K, compared with 11,400K from the spectra. 
More observations of the star WDJ 0303-0808 are required, 
both with high SNR spectroscopy and time-resolved pho- 
tometry, for a better understanding of its structure. 

5. Concluding remarks 

We report the discovery of seven new pulsating DA stars, 
bringing the total to 126 known ZZ Cetis in the narrow 
temperature range 12 270 K > T e s > 10 850 K, correspond- 
ing to the existence of a partial ionization of hydrogen and 
development of a sub-surface convection zone. As two of 
these stars were previously reported as NOV, we need to 
lower the detection limits for all other NOVs observed by 
Mukadam et al. (2004) and Mullally et al. (2005), map the 
instability strip in both T e ff and log g using accurate de- 
terminations of these parameters for the new variables and 
NOVs, for which wc need higher SNR than those achieved 
in the SDSS spectra. 
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